Decoherence from Spin Environments 
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We examine two exactly solvable models of decoherence - a central spin-system, (i) with and (ii) 
without a self-Hamiltonian, interacting with a collection of environment spins. In the absence of a 
self-Hamiltonian we show that in this model (introduced some time ago to illustrate environment- 
induced superselection) generic assumptions about the coupling strengths can lead to a universal 
(Gaussian) suppression of coherence between pointer states. On the other hand, we show that when 
the dynamics of the central spin is dominant a different regime emerges, which is characterized by 
a non-Gaussian decay and a dramatically different set of pointer states. We explore the regimes of 
validity of the Gaussian-decay and discuss its relation to the spectral features of the environment 
and to the Loschmidt echo (or fidelity). 

PACS numbers: 03.65.Yz;03.67.-a 



I. INTRODUCTION 

A central spin-system S interacting with an environ- 
ment £ formed by N independent spins through the 
Hamiltonian 
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N 

fc=i 



(1) 



may be the simplest solvable model of decoherence (we 
use the standard notation according to which (7^^^^ and 
(Ti, i = x,y,z, denote Pauli operators acting on the k- 
th environmental spin and on the central system) . This 
Hamiltonian was studied some time ago |1| as a simple 
model of decoherence. It was used to show that rela- 
tively straightforward assumptions about the dynamics 
can lead to the emergence of a preferred set of pointer 
states due to environment-induced superselection (eins- 
election) 0, |2| ■ Such models have gained additional im- 
portance in the past decade because of their relevance to 
quantum information processing 3] . 

The model described by |^ was particularly useful to 
ilustrate the nature of decoherence in the context of a 
measurement. In such case the central spin is used as a 
simple (two state - one bit) approximation for a mem- 
ory of classical apparatus. Then, it is natural to neglect 
the effect of the system's self-Hamiltonian. As a conse- 
quence, the eigenstates of the interaction Hamiltonian 
emerge as preferred pointer states of the system (defined 
as the ones which are "least perturbed" by the interac- 
tion with the environment 2]). Thus, the eigenstates of 
the az operator (denoted here as |0) and with eigen- 
values and —1 respectively) are dynamically selected 
by the interaction with the environment. Indeed, these 
states are not perturbed by the interaction while other 
superpositions rapidly decay into their mixtures. 

Neglecting the self-Hamiltonian of the system is not al- 
ways a reasonable approximation. Studies of decoherence 
without such assumption have also been carried out us- 
ing, mostly, the Quantum Brownian Motion as a paradig- 



matic example In such case, pointer states do not 
coincide with the eigenstates of the interaction Hamilto- 
nian but can range from coherent states for the QBM case 
to eigenstates of the system's Hamiltonian iQ] . Their 
properties are determined by the interplay between the 
self-Hamiltonian and the interaction with the environ- 
ment. 

In this paper we will study a generalization of the 
above simple model described by the Hamiltonian 



Ht = Act, 
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(2) 



This simple model includes both the effect of the evolu- 
tion of the central system and its coupling with the spin 
environment. 

The purpose of our study is twofold. First, in Section 
II we will analyze again the case where the central spin 
has no self-Hamiltonian (A = above). Our goal is to 
show that - with a few additional natural and simple as- 
sumptions about the distribution of coupling strengths 
gk - one can evaluate the exact time dependence of the 
reduced density matrix of the central spin. In fact, we 
will demonstrate that the off-diagonal components dis- 
play a Gaussian (rather than exponential) decay. In this 
way we will exhibit a simple soluble example of a situ- 
ation where the usual Markovian R| assumptions about 
the evolution of a quantum open system are not satisfied 
at any time. 

Then, in Section HI we will consider the more complex 
case with non-trivial dynamics (A 7^ above). We will 
show that, under the same natural assumptions made in 
Section II about the distribution of coupling strengths 
in the interaction Hamiltonian, the problem can also be 
solved exactly. The solution will enable us to study two 
very important features of the decoherence process. We 
will analyze the nature of pointer states and also the 
way in which the reduced density matrix of the central 
spin evolves in time. In this case, the decay of the off- 
diagonal component is not Gaussian but displays long 
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time algebraic tails which we obtain analytically. The 
most probable pointer states will be shown to range from 
eigenstates of cr^ in the small A limit to eigenstates of 
ax in the opposite limit of large A (result that can be 
expected based on the considerations presented in |^). 
In Section IV we will summarize our results which, apart 
from their implications for decoherence, could also be 
relevant to quantum error correction [8| where precise 
knowledge of decoherence is essential to select an efficient 
strategy to defeat it. 



II. STATIC SYSTEM - GAUSSIAN 
DECOHERENCE 

Here we will consider the system described by Eq. Q . 
We begin by outlining how to solve this model exactly, 
and how to find the time dependence of the elements of 
the reduced density matrix of the system. Let us consider 
an initial state for the combined system-environment of 
the form 



|*5£(0)) = (a|0) 



2™-l 

n=0 



(3) 



Here \n) are the states of the computational basis of the 
environment that diagonalizes Tiss . The k-th digit of the 
binary form of n, n^, represents the state up or down in 
the Z axis of the fc-th spin of the environment. The main 
assumptions above are that the initial state is a product 
(no initial entanglement between the system and environ- 
ment) and that the total state is pure. Both conditions 
can be easily relaxed, but choosing Eq. (|3Jl simplifies the 
presentation. The state of S£ at an arbitrary time is 
given by: 

\^seit))^a\0)\£o{t))+b\l)\£^{t)), (4) 



with 



2~-l 

= |fi(-i)>, (5) 



and where 



N 

Bn^Y.^-ir9k (6) 
The reduced density matrix of the system is then: 

PS = Tr£|*5£(t))(*5£(i)| 

= |anO)(0|+a6V(t)|0)(l| 

+ a*for*(i)|l)(0| + |feni)(l|, (7) 

where the decoherence factor r{t) ~ {£i{t)\£o{t)) can be 
readily obtained: 



it) - E 



(8) 



It was shown in (using some simplifications to be 
discussed below) that for t > 0, r{t) decays rapidly to 
zero, so that the typical fluctuations of the off-diagonal 
terms of ps will be small for large environments. There- 
fore, the decoherence factor tends to zero (Ir(i)P') — > 

0, leaving ps approximately diagonal in a mixture of the 
pointer states {|0) , |1)} which retain preexisting classical 
correlations. 

We will show in this section that, for a fairly generic set 
of assumptions, the form of r{t) can be further evaluated 
and that - quite universally - it turns out to be approx- 
imately Gaussian in time. To prove this we will only 
require that the couplings gt of Eq. are sufficiently 
concentrated near their average value so that their stan- 
dard deviation {^{g^ — {gkiY') exists and is finite. When 
this condition is not fulfilled other sorts of time depen- 
dence become possible. In particular, r{t) may be expo- 
nential when the distribution of couplings is, for example, 
Lorentzian. 

To obtain our result we rewrite Eq. (|SJ) as 



r{t) 



-iBt 



r/{B)dB, 



(9) 



that is, the decoherence factor is the Fourier transform 
of a characteristic function 



2"-l 

r]{B) = E \CnfSiB-Br,). 

n=0 



(10) 



Eq. (|10|l is a particular case of the more general strength 
fimction or local density of states 0, 

,7(S)-El('^>l*5^)l''^(^-^>)' (11) 



where are the eigenfunctions of the full Hamiltonian 
with eigenenergies Bx. 

The discussion of decoherence in our model is thus di- 
rectly related to the study of the characteristic function of 
the distribution of couphng energies ri{B). Since the B^s 
are sums of gk^s (that we assume independent of each 
other), equation @ makes r{t) itself a product of charac- 
teristic functions of the distributions of the couplings gt- 
Thus, the distribution of _B„ belongs to the class of the 
so-called infinitely divisible distributions [T]| . The 
behavior of the decoherence factor r{t) - characteristic 
function of an infinitely divisible distribution - depends 
only on the average and variance of the distributions of 
couplings weighted by the initial state of the environment 
EkiiJ- 

Assuming that the variance of the couplings gk is fi- 
nite, we claim that for reasonable assumptions on the 
initial state of the environment (the coefficients c„), and 
N sufficiently large, r]{B) has in general a Gaussian form. 
Therefore, the decoherence factor decays as a Gaussian 
with time. We will show this behavior with some exam- 
ples where an exact solution is possible, and discuss the 
regime of validity of the conjecture. 



Let us consider first the simplest case where all cou- 
plings are equal, = g, and all the spins of the environ- 
ment have the same initial state, 

N 

1*5^(0)) = (a |0) + 6 |1)) (g)K |0), + (12) 

k=l 

with ak — a and (3k — P for all k. The decoherence factor 
then takes the simple form r(t) = (|ape*f*-f |/3|2e~^9*)^. 
Expanding this expression we obtain 

^(*) = E ( fc) lal''!/?!'^"^''^' exp Hgt{2k - N)]. 

_(13) 

As follows from the Laplace-dc Moivrc theorem [l^l, for 
sufhciently large the coefficients of the binomial ex- 
pansion of Eq. (|13|l can be approximated by a Gaussian, 
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(14) 



Therefore for large 



exp 



77(i?) 



8Ar|a/3|2 



v/87rAf|a/3|2 



(15) 



This generic behavior can be interpreted as a result of the 
law of large numbers JiJ]: the energies -B„ of the com- 
posite S£ system can be thought of as being the termi- 
nal points of an A^-step random walk. The contribution 
of the fc-th spin of the environment to the random en- 
ergy is +g or —g with probability |ap or respectively 
[Fig. (^a)]. Therefore, the set of all the resulting energies 
must have an (approximately) Gaussian distribution. 

We can carry out the same argument in the more gen- 
eral case of Eq. H12|l for different couplings and initial 
states for the spins of the environment. Here, 
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r{t) = Y[{ 



+ \Pk 



-iQkt 



)• 



(16) 
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The "random walk" picture that yielded the distribu- 
tion of the couplings remains valid [see Fig. (^b)]. How- 
ever, now the individual steps in the random walk are no 
longer all equal. Rather, they are given by the set {gk} 
[see Eq. (Q] with each step gk taken just once in a given 
walk. There are 2-^^ such distinct random walks Wn, one 
for every state \n) of the environment. Each walk con- 
tributes to r]{B) with the weight given by the product of 
the relevant jafcj^ and |/3/c|^, or right (fc € W^+) and left 
{k G W~) "steps" respectively. The weight of the n-th 
walk is then given by 



(17) 





FIG. 1: The distribution of the energies obtains from the ran- 
dom walks with the steps given by the coupling size and in 
the direction {+gk or —Qk) biased by the probabilities jofej^ 
and |/3fcp as in Eq. IIUII (although in these examples we set 
Qij^ — 1/2). (a) When all the couplings have the same size 
Qk = g (Eq. 1131 '). a simple Newton's triangle leads to an ap- 
proximate Gaussian for the distribution of energies, (b) When 
the couplings differ from step to step (Eq. I21ll 'l. the resulting 
distribution still has a approximately Gaussian envelope for 
large A^. 



The terminal points Bn of the random walks may or may 
not be degenerate: As seen in Fig. (Q), in the degenerate 
case, the whole collection of 2^ random walks "collapses" 
into A -|- 1 terminal energies. More typically, in the non- 
degenerate case [also displayed in Fig. (Q)] , there are 2^ 
different terminal energies In any case, the "enve- 
lope" of the distribution r]{B) will be Gaussian, as we 
shall argue below. 

Let us compute the characteristic function ri{B). If we 
denote Xk the random variable that takes the value +gk 
or —gk with probability jafcj^ or respectively, then 
its mean value ak and its variance bk are 



Ofe = {\akf - \Pkf)9k, 
bl = gl-al=A\ak\^\f3k 



I'd- 



(18) 



The behavior of the sums of A^ random variables Xk (and 
thus, of their characteristic function) depends on whether 
the so-called Lindeberg condition holds j^. It is ex- 
pressed in terms of the cumulative variances sff — 
and it is satisfied when the probability of the large indi- 
vidual steps is small; e.g.: 



P{ max \gk - a^j > tsn) 

l<k<N 



N- 



0, 



(19) 



for any positive constant r. In effect, Lindeberg condition 
demands that sn be finite: when it is met, the resulting 
distribution of energies B — ^Xk is Gaussian 



B- B 
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sn 



< X 



N- 



dy, (20) 
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FIG. 2: Assumed distribution of the couplings gt, and re- 
sulting distribution of the eigenenergies -B„ (left panels) for 
N = 6 and = 24. In the case of \akf = 1/2 this distri- 
bution is in effect the "strength function" (local density of 
states). The corresponding decoherence factors r{t) for dif- 
ferent initial conditions with A = 6 (dashed lines), A'' = 24 
(thin solid lines) and the average (bold line) is shown on the 
right. Note the rapid convergence to a Gaussian behavior. 



where Bn — J2k'^k- This imphes 



: exp 



2s2 



(21) 



an expression in excellent agreement with numerical re- 
sults already for modest values of N. After applying 
the Fourier transform of Eq. @, this distribution of 
energies yields a corresponding approximately Gaussian 
time-dependence of r{t) [Fig. 



r{t) 



(22) 



Moreover, at least for short times of interest for, say, 
quantum error correction, r(t) is approximately Gaussian 
already for relatively small values of N. This conclussion 
holds whenever the initial distribution of the couplings 
has a finite variance. Note that in particular, we did not 
have to assume "randomness" of the couplings gk [see 
e.g. Eqs. inj-lini)]. 

A random initial state for the environment (not neces- 
sarily a product state) instead of Eq. H12|) gives basically 
the same result. In this case typically c„ ~ 2~^/^e^'^" , 
with (/)„ a random phase between and 27r. From Eq. 

(Cnj, 
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fe=0 ^ ^ 

and, as above, the Gaussian limit for large N applies. 

It is also interesting to investigate cases when Linde- 
berg condition is not met. Here, the possible limit distri- 
butions are given by the stable (or Levy) laws jljj . One 
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FIG. 3: Same as Fig. 2 but for a Lorentzian distribution of 
the couplings gt- In this case r{t) decays exponentially. The 
histogram and the dashed line in r(t) correspond to A" = 20, 
the straight thin line is a particular case for A" — 100 and 
the thick line is the average. We note that the convergence 
is slower than in the Gaussian case of Fig. 2, because real- 
izations of gk are more likely to have one or two dominant 
couplings. Therefore, although the average shows a clear ex- 
ponential decay, fluctuations are noticeable even for large A". 



interesting case is a Lorentzian distribution of couplings, 
which yields an exponential decay of the decoherence fac- 
tor [see Fig. Q]. Such a distribution could be obtained 
for instance by considering dipolar interaction between 
spins randomly placed in a sample. The long range na- 
ture of the interaction gives rise to the Lorentzian distri- 
bution and therefore to the exponential decay that can 
be deduced by statistical arguments 



A. Relation to the Loschmidt echo 

The Fourier transform of the strength function ri{B) is 
also related to the Loschmidt echo 13] (or fidelity) in the 
so called Fermi Golden rule regime. The fact that the pu- 
rity and the fidelity have closely related decay rates has 
been recently shown for the case of a bath composed 
of non-interacting harmonic oscillators. In this sense our 
results could be interpreted as an extension of the discus- 
sion of Ref. ^3 to spin environments. 

The connection with fidelity is more easily seen if we 
write a generalized version of the Hamiltonian (^), 



(|0) (0|®H° 



ID (II 



(24) 



The decoherence factor is then the overlap of the initial 
state of the environment [^'^(O)) evolved with two differ- 
ent Hamiltonians, 

r(t) = (*£(0)| e*'"?*/2g-mi t/2 \-qj^^Q)) ^ (25) 

which clearly has the form of the amplitude of the 
Loschmidt echo for the environment with the two states 
of the system as the perturbation. In the model of Eq. 
O, T^g = -nl and thus 



r(t) 



(vI/£(0)|e-^«^*|vI/£(0)) 



(26) 



This expression is the survival probability of the initial 
state of the environment under the action of the Hamil- 
tonian Ti.^, which is known to be the Fourier transform 
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of the strength function |^^. This connection provides 
another way to understand Eq. . 



III. DECOHERENCE AND DYNAMICS 

In this section we will study the more general Hamilto- 
nian of Eq. ©, that is we will include a self Hamiltonian 
to the central system. The results of the previous sec- 
tion will be contained in the limit of A = 0, however we 
will see that for any finite A the behavior of the decoher- 
ence factor will be non-trivially different from what we 
obtained in the previous section. 

Despite its more complex appearence, the model given 
by Eq. ||2Jl is still exactly solvable Since the states 
\n) {n\ of the environment commute with the Hamilto- 
nian, we can write the evolution operator for the com- 
bined system-environment as 



(27) 



n=0 



with 



UB„{t) = /cos(ri„t) _ , (^^-^" + ^"^) sin(r!„t), (28) 



and rj^j = A^ -I- _B^. The physical interpretation of this 
results is that for every state of the environment \n) the 
effective dynamics of the system is given by a magnetic 
field fin = (A, 0, Bn) in the XZ plane. Seen from this 
perspective, the decoherence is produced by the disper- 
sion of the fields i?„. 

The reduced density matrix of the system at an arbi- 
trary time t is 

2"-l 



P{t)= E \^n\^UBAt)p{mlSt). (29) 

ri=0 

or, transforming the notation and using Eq. UlUf) . 



p{t) - / UB{t)p{0)ul{t) rj{B) dB. 



(30) 



For simplicity, we will work with the polarization vec- 
tor p, such that p = {I + p - ct)/2. Thus, 



p(t,B)fi{B) dB 



For an arbitrary time t, we find 



(31) 



, ,A'^ + B'^cos(2nBt) , , B , ^ , , ,2AB , 

p,{t, B) = p,(0) „2 ^ ' - Py{<d)—sin{2nBt) +P^(0)— ^ sin2(r!Bt), 

\Ib ^^B 



(32a) 



Pyit, B) = py{0) cosi2nBt) + [p^(o)B - ApM] 



(32b) 



p,{t,B)^pM 



B^ + A2 cos(2rJBt) 



^2 



2 A_B 

(0) 7]2- sin' inBt)+ py (0) ^ sin(2r!si) , 



A 

'^B 



(32c) 



According to the results of the previous section, in gen- 
eral for large TV we can assume a Gaussian shape for ri{B) . 
By using a Gaussian centered around zero. 



1 



exp {-B'^/2t 



N ' 



(33) 



Eqs. (|32|l simplify because the odd terms in B don't 
contribute to the final result. 

Using these assumptions, we were not able to obtain a 
solution of the integral in Eq. (|31|l for arbitrary values of 
SAT and A. However, we can solve the two limiting cases 
SAf ^ A and sat <C A, which turn out to give non-trivial 
results. 

Let us consider first the case where sat ^ A, that 
is, where the central spin dynamics is so slow that its 



behavior should approach that obtained in the previous 
section. Indeed, for short times [t <^ A^^), using a Tay- 
lor expansion of Eqs. (|32|l around A = one finds 

P.{t) - P.(0)e-^*'^" 

Py{t) = p^(0)e-2*'-« _p^(0)A /jErf(\/2sArt) 

sn s ^ 

Pz{t) ^ p,iO)+Py{0) — J^ETf{V2sNt), (34) 
Sat V ^ 

where Erf(a;) is the error function. To obtain the long 
time behavior, we need to perform the integrals on B by 
stationary phase approximation. In the limit t ^ A~^ 
we find 
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FIG. 4: Behavior of the components of the polarization vector 
for long times for A/sn = 0.1 (top panel) and A/sn = 5 
(bottom panel) . Numerical results in solid lines and analytical 
predictions in dashed lines. 



Px{t) 
Py{t) 
Pz{t) 



— + , cos {2At+ — 



Py(0) COS (2 At + -PziO) sin (2 At + ^ 



P.(0) 



1-7 



A 



V2. 



sn 



-4- cos f 2Ai + ^ 
2s%t V 4 



sin f 2At + ^ 

2sjft V 4 



(35) 



with 7(x) = y^xe^ (1 — Erf(x)). In this limit, 

Note that for any A ^ the X component of the 
polarization does not decay to zero, indicating that the 
decoherence process is not completely effective in this 
direction [the Y component does go to zero for large times 
due to the symmetry of Hamiltonian Q]. Also, note that 



even a small self-Hamiltonian of the system always ends 
up turning a fast (Gaussian) decay into a slow (power 
law) one. 

In the opposite limit of strong self-dynamics of the sys- 
tem. Sat ^ A, we can obtain an expression valid for all 
times by expanding Hb ~ A -f- B^/2A. After some alge- 
bra. 



Pz{t) - pAO) 



1/4 



Py (0) COS ( 2 At + - arctan ^ 



Pz (0) sin ( 2 Ai + ^ arctan 



SN 



COS 2At - 



i arctan 



(l + ^) 



sin ( 2Ai 



i arctan 



(l + ^) 



1/4 



(36) 




FIG. 5: Schematic solution of the Hamiltonian ((SJ for two 
initial polarizations in the Z axis (left panel) and in the X 
axis (right panel) . Supposing that B takes only two possible 
values, ±S]v, the solution for each field B is the precession of 
the polarization around il± = {A,0,±sn)- The total polar- 
ization is the average of the two cones, which gives a small 
residue along the Z axis. The polarization in X is almost 
fully preserved. 



In the long time limit (t 3> A/s^) these expressions are 

equal to Eqs. (|35|l . only that now 7 ^ ^ ~ 1. The 

results above for large and small A agree well with nu- 
merical simulations, as shown in Fig. Q). 



A. Pointer basis 

The above results allow us to draw some conclusions 
about the nature of the decoherence process and the 
pointer states which are dynamically selected by the en- 
vironment. First, we can notice that for long times the 
polarization vector converges to a certain value (which, 
in general, depends on A and other parameters of the 
model). Second, we note that when A ^ the X com- 
ponent of the polarization vector does not decay to zero 
but is resilient to the interaction with the environment. 
This is the case even if the system interacts with the 
environment through the Z component of the spin. 

The states which are dynamically selected by the en- 
vironment are dramatically different in the two oposite 



regimes we examined above. For small values of A, the 
eigenstates of the Z component of the central spin are 
pointer states. They are minimally perturbed by the in- 
teraction with the environment (in the previous section, 
where A = was assumed, this emerged as an exact re- 
sult since Pz is conserved). However, for large A (i.e. 
A ^ sn), the fact that ^^^(i ^ cxd) ~ 1 is a signature of 
the decoherence process selecting a completely different 
set of pointer states. In fact, in this case, the pointer 
states turn out to be eigenstates of the system Hamil- 
tonian, which is proportional to a^. Thus, this model 
enables us to examine these two very different situations: 
one where the interaction with the environment domi- 
nates (A <C Sn) and Cz eignestates are selected; the 
other where the self-Hamiltonian of the system domi- 
nates (A ^ Sn) and ax eigenstates are selected. 

The regime where the Hamiltonian of the system domi- 
nates (or, more precisely, where the environment is much 
slower than the system) was analyzed in a more gen- 
eral contex before and has a natural interpretation 
here: This behavior is simply the one corresponding to 
the strong decoupling regime observed in Nuclear Mag- 
netic Resonance 0|. There, the presence of a strong 
magnetic field in the XorY axes effectively decouples the 
spectrally resolvable spins of a sample (whose interaction 
is Z dominant). The standard picture of this decoupling 
regime is that by rotating the polarization rapidly enough 
around X, any interaction in another axis is strongly sup- 
pressed and the spins effectively "decouple" . 

There is an instructive physical picture to understand 
these results. Instead of using a continuos distribution for 
B, let us suppose that B can only take two values, ri{B) — 
[d{B - Sn) + S{B + sn)]/2, with sn < A. The classical 
solution for the evolution of the polarization vector is 
the precession of p around ft± — (A, 0,±SAr), as shown 
in Fig. (0) for two possible initial conditions of p. The 
polarization vector of the reduced system is the average 
of the two cones corresponding to the precession around 
r2+ and n_ The presence of a small component A in the 
Z axis tilts the precession cones so that their average is 
almost 1 in the X direction and has a small residue on the 
Z axis (the Y component cancels due to the symmetry). 
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FIG. 6: Bloch sphere representation of our main results. Pos- 
sible initial states at t = are represented on the left by a full 
sphere. Intermediate times (center) are similar for A = and 
A ^ Sat , where decoherence reduces the Bloch sphere equally 
fcist (Gaussian) in the X and Y axes. For large A, Y and Z 
axes are decohered in a slower algebraic way. For long times, 
right panel, the A = case is completely decohered to the Z 
axis, while the small A retains some polarization along the 
X axis. The large A case is almost the opposite, retaining 
almost all polarization in the X axis and a small residue in 
the Z axis. 

IV. CONCLUSIONS 

We have studied a very simple model of decoherence 
due to spin environments. We showed that the decoher- 
ence factor will generically have a Gaussian decay when 
there is no self-Hamiltonian for the system. We note that 
similar behavior was observed for short times by studying 
the decoherence process in models where the largest en- 
ergy scale is the system-environment interaction strength 
^M- A model similar to ^ is used in the NMR setting 
[l^l to compute corrections to the second and fourth mo- 
ments of the decay of the polarization signal. Here the 
idea is to treat the interaction with surrounding spins 
as an effective local magnetic field that shifts the Lar- 
mor frequency inhomogenuosly across the sample. The 
statistical treatment used in contrasts with the ex- 
act solution presented in this work, even in the pres- 
ence of a self-Hamiltonian of the central spin. Thus, our 



model has applicability and relevance to a larger class of 
physical situations. There is a substantial body of work 
[ITL I20I I21L l2a | on decoherence due to spin environments, 
stimulated in part by the interests of quantum compu- 
tation. Our results are most relevant for quantum error 
correction and other strategies to fight decoherence in a 
quantum computer. For example Gaussian time depen- 
dence of the decoherence factor would suggest a different 
(more frequent) error correction than the exponential de- 
pendence often assumed with little or no justification. 

We also showed how by adding a self-Hamiltonian for 
the system one can dramatically change the main features 
of the decoherence process. Even for the case of slow dy- 
namics of the system, we found that for long times the 
initial Gaussian behavior changes to a power law. On the 
other end, when the self-Hamiltonian is much stronger 
than the interaction with the environment, the whole 
process changes its nature. The decay is predominantly a 
power law. Moreover, the pointer states now correspond 
to eigenstates of the system rather than eigenstates of 
the system-environment interaction [^. For illustrative 
purposes, our results are summarized schematically in 
Fig.® using the Bloch sphere representation. 

Our results, though interesting, arise from a very sim- 
plified model. A logical step for future research is the 
inclusion of intra-bath interactions. The entanglement 
thus created between spin baths will surely have an im- 
pact on the amount of decoherence in the system 23] . 

Possible experimental applications of our considera- 
tions are in nuclear magnetic resonance, and in any other 
situation where two-level systems interact with spin envi- 
ronments. Another area of impact of our results is in the 
characterization of the process that leads to redundancy 
in the environment of the classical information about the 
system [2j| . The relation between the decoherence factor 
and the strength function might prove useful in the phys- 
ical setting of strongly interacting fermions, where it has 
been shown that the strength function takes a Gaussian 
shape . It is our hope that the simple analytic model 
described here will assist in gaining further insights into 
these fascinating problems. 
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